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If dark energy (DE) couples to neutrinos, then there may be apparent violations of Lorentz/ CPT 
invariance in neutrino oscillations. The DE-induced Lorentz/ CPT violation takes a specihc form 
that introduces neutrino oscillations that are energy independent, differ for particles and antipar- 
ticles, and can lead to novel effects for neutrinos propagating through matter. We show that 
ultra-high-energy neutrinos may provide one avenue to seek this type of Lorentz/ CPT violation in 
Uf_i-VT oscillations, improving the current sensitivity to such effects by seven orders of magnitude. 
Lorentz/ CPT violation in electron-neutrino oscillations may be probed with the zenith-angle depen- 
dence for high-energy atmospheric neutrinos. The "smoking gun," for DE-neutrino coupling would, 
however, be a dependence of neutrino oscillations on the direction of the neutrino momentum rel- 
ative to our peculiar velocity with respect to the CMB rest frame. While the amplitude of this 
directional dependence is expected to be small, it may nevertheless be worth seeking in current data 
and may be a target for future neutrino experiments. 
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I. INTRODUCTION 



The accelerated cosmic expansion \U poses difficult 
questions for theoretical physics [1, H, Is it simply 
due to a cosmological constant? Is some new negative- 
pressure dark energy (DE) required? Is general relativ- 
ity modified at large distance scales? The major thrust 
of the empirical assault on these questions has been to 
determine whether the expansion history and growth of 
large-scale structure are consistent with a cosmological 
constant or require something more exotic [H]- 

However, it may be profitable to explore whether there 
are other experimental consequences of the new physics — 
which we collectively refer to as DE, although it may in- 
volve a modification of gravity rather than the introduc- 
tion of some new substance — responsible for accelerated 
expansion. If cosmic acceleration is due to a cosmolog- 
ical constant (i.e., if general relativity is valid and the 
equation-of-state parameter \s w = —1), then the vac- 
uum is Lorentz invariant. If, however, something else is 
going on, then the "vacuum" has a preferred frame: the 
rest frame of the cosmic microwave background (CMB). 
If, moreover, dark energy couples somehow to standard- 
model particles, then there may be testable (apparent) 
violations of Lorentz invariance. For example, if DE is 
coupled to the pseudoscalar FF of electromagnetism ^ , 
there may be a "cosmological birefringence" that rotates 
the linear polarization of cosmological photons; CMB 
searches for such a rotation ^ constrain this rotation 
to be less than a few degrees 8] . 

Here we explore DE-induced Lorentz/ CPT- violating 
effects in the neutrino sector. We show that the form 
of a Lorentz-violating coupling between neutrinos and 
dark energy is highly restricted under fairly general as- 



sumptions.^ The coupling engenders an additional source 
for neutrino mixing (e.g., Ref. (To|). resulting in neu- 
trino oscillations with a different energy dependence 
than vacuum oscillations and different oscillation prob- 
abilities for neutrinos and antineutrinos. While similar 
Lorentz / CPT- violating oscillations have been considered 
before PEIll, we emphasize here that cosmic accel- 
eration dictates a specific form for such effects. 

Data from Super-Kamiokande and K2K 14] and 
AMANDA/IceCube [15] already tightly constrain CPT- 
violating parameters for v^-v^ mixing, and those from 
solar-neutrino experiments and KamLAND ^16,] do so 
for v^-v^j^ mixing. However, the effects of DE-induced 
CPT violation become more significant at higher energies 
[l3]. Here we show that next-generation measurements 
of ultra-high-energy neutrinos produced by spallation of 
ultra-high-energy cosmic rays will increase the sensitiv- 
ity to CPT-violating v^-Vt oscillations by seven orders 
of magnitude. We also show that these CPT-violating 
couplings may lead to novel effects in the zenith-angle 
dependence for atmospheric neutrinos in the ^100 GeV 
range. 

While such CPT-violating effects, if detected, could be 
attributed simply to intrinsic CPT violation in funda- 
mental physics, not related to DE, a DE-neutrino cou- 
pling further predicts a directional effect: the neutrino- 
mixing parameters depend on the neutrino propagation 
direction relative to our peculiar velocity with respect 
to the CMB rest frame. While this signature will likely 
remain elusive even to next-generation experiments, it 
would, if detected, be a "smoking gun" for DE beyond a 



The coupling of neutrinos to dark energy has also been considered 
in the context of "mass- varying neutrinos" [^, but that imple- 
mentation of the DE-neutrino coupling does not lead to the type 
of Lorentz/ CPT-violating effects we discuss here. 
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cosmological constant. It is therefore worth considering 
as a long-range target for future neutrino experiments. It 
may also be worthwhile to search current data in case an 
implementation of DE-neutrino coupling different from 
that we consider here leads to a different energy depen- 
dence for these directional effects. We therefore work 
out explicitly the directional dependence to aid experi- 
mentalists who may wish to look for such correlations in 
current data. 

Below, we first derive in Sec. |TI] the form of 
the Lorentz/CPT violation allowed by a DE-neutrino 
coupling and discuss the resulting neutrino-oscillation 
physics. In Sec. IIIII we apply the formalism to cosmo- 
genic ultra-high-energy neutrinos, and obtain projected 
sensitivities of future detectors to these effects in Vf^-Ur 



oscillations. In Sec. llVl we discuss matter- induced effects 
for Vf, oscillations in high-energy atmospheric neutrinos 
in the presence of Lorentz-invariance- violating mixings. 
Concrete formulas for the directional dependence on os- 
cillation probabilities are given in Sec. |Vl Finally, we 
discuss some theoretical implications in Sec. lVII and sum- 
marize and conclude in Sec. IVIII 



II. THE DARK-ENERGY-NEUTRINO 
COUPLING 

A. General Formalism 

Following Ref. fisj, the neutrino fields are denoted 



by Dirac spinors {vei^ 
jugates by {; 



} and their charge con- 
}, where z/^^c = v'^ = 
Cv^ is the charge-conjugated spinor, and C is the 
charge-conjugation matrix. The 2N fields (where N 
is the number of flavors) and their conjugates are ar- 
ranged in a single object va where A ranges over 
e,M,r, ••• ,e'=,^'=,r'=, • • • . 

With a canonical kinetic term in the neutrino La- 
grangian, the most general Lorentz/ CPT- violating Dirac 
equation is,^ 



{l^^d^ ~ MAB)yB = 0, 



(1) 



where 



Pr = (l+75)/2. The 27V X 2iV" mass matrix m^j is written 
in terms oi N x N matrices L, R, and D, through 

Here, R and L are the right- and left-handed Majorana 
neutrino masses (L = is required if electroweak gauge 
invariance is preserved), and D is the Dirac-mass matrix. 
The R and L matrices are required to be symmetric, and 
R, L, and D can most generally be complex. 



B. Dark-energy— induced Lorentz violation 

Lorentz violation in Eq. ^ is parametrized by the 
four- vectors a^, b^, and the antisymmetric tensor i/^". 
The parameters and 6^ are both CPT and Lorentz vi- 
olating, while H^^'^ is Lorentz violating but CPT conserv- 
ing. While these parameters are non-zero for the most 
general Lorentz/ CPT- violating Dirac equation [l^, the 
allowable forms for a^, b^, and H^'^ are highly restricted 
if the Lorentz/ CPT violation is induced by coupling to 
dark energy. 

The smallness of the CMB quadrupole demands that 
the 3-dimensional hypcrsurfaccs of constant DE density 
must be closely aligned with those of constant CMB tem- 
perature [Tsl. The preferred frame associated with the 
cosmic expansion is then parametrized by a unit four- 
vector which is orthogonal to surfaces of constant 
CMB temperature; i.e., in the CMB rest frame, it is 

= (1, 0, 0, 0). The symmetry of the problem thus dic- 
tates that (X Z'" and b^' oc Z'". The tensor H^'' is 
antisymmetric, and there is no way to construct an anti- 
symmetric tensor H^'^ from a single four- vector; we thus 
expect H^"^ = for DE-neutrino coupling. 

Furthermore, since neutrinos are produced and in- 
teract in weak eigenstates, it is only the combination 
(^i)ab = ('^ + b)ab that is relevant for neutrino phe- 
nomenology. Thus, the Lorentz/CPT violation induced in 
neutrino physics can be parametrized entirely by a single 
four-vector-valued (aL)^}, oc l^^ matrix in the flavor space. 



M 



AB 



mAB 
1 



im^ABlb 



-H 



b^. 



(2) 



Neutrino Oscillations 



The usual mass terms are m -\- im^j^ = m^PL -\- mnPn, 
where mji = (tol)^ — m -\- im^, Pl = {1 — 75)/2, and 

^ Additional possibilities arise with a non-canonical kinetic term; 
we comment briefly on possible consequences below. 



The propagation of the flavor eigenstates is then de- 
scribed by an effective Hamiltonian 



J 



{heS, 



ab 



_ ( pSab + {rn^)ab/'2p+ {aL)abPtJ./P 





pSab + {rh^)lJ2p - {aL)llPt,/p 



(4) 
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where the flavor indices a and b run over the flavor eigen- 
states e, /j,, and r and e"^, /i'^, and t'^. Here, p = with 
p the neutrino momentum, and rh^ = mimj is the usual 
mass matrix, with mi = L — DR~^D^ . 

Equation (|4]) has several implications: (i) Since the 
matrix is block-diagonal, there is no mixing between neu- 
trinos and antineutrinos (as may arise in more general 
Lorentz-violating scenarios [isj). (ii) Since ql appears 
with opposite sign in the neutrino and antineutrino en- 
tries in the Hamiltonian, a nonzero implies (apparent) 
CPT violation — i.e., the propagation of neutrinos and an- 
tineutrinos is not the same. Thus, for example, if the 
anomalous LSND results had stood, the CPT-violating 
explanations (e.g., Ref. jT^) for them [13] may have im- 
plied DE-neutrino coupling, (iii) The mixing induced by 
DE-neutrino coupling is energy independent (like in the 
Mikheev-Smirnov-Wolfenstein, or MSW, effect [2ll|). as 
opposed to vacuum mixing, which declines as E~^. Thus, 
these effects will become increasingly visible at higher 
energies. The detailed form of CPT violation implied 
by this effect is also thus different than that obtained 
with different Am^ for neutrinos and antineutrinos. (iv) 
There may also be novel effects for neutrinos propagating 
through matter, an effect we discuss further in Sec. IIVI 
below. 



Finally, (v) the neutrino oscillations induced by DE- 
neutrino coupling are frame dependent. If the observer 
is in the rest frame of the CMB, then [a^Yp^ cx E, and 
neutrino oscillations are independent of the neutrino di- 
rection. However, the Solar System moves with respect 
to the CMB rest frame with a velocity v ~ 370 km s^^. 
DE-induced neutrino oscillations will therefore depend 
on {aL)^p^ cx E{\ — V ■ p)^ where p is the neutrino- 
propagation direction and v is our peculiar velocity with 
respect to the CMB rest frame. There will thus be an an- 
nual modulation in solar-neutrino oscillations, a diurnal 
modulation in laboratory neutrino- mixing experiments, 
and a direction dependence in oscillations of cosmogenic 
neutrinos. 

Since neutrino mixing arises only as a consequence of 
the traceless part of the propagation Hamiltonian, the 
DE-neutrino coupling must (like the vacuum mass ma- 
trix) be flavor-violating if neutrino oscillations arc to be 
affected. 



D. Two- flavor oscillations 

The evolution equation for DE-induced two-flavor mix- 
ing is of the form. 



dt 



1 
2 



cos 2^1, - mcfr(l - v ■ p) cos26d 
^ sin 29., + mcff(l - V • p) sin 261^6-"' 



Am' 
2E 
Am 
2E 



sin 2f?„ - 
COS 29jj 



- moff(l — V ■ p) sin 20^6"' 
+ TOcff (1 ~ V ■ p) COS 29d 



(5) 



where rucs is an effective mass parameter, and 9d and rj 
are a mixing angle and relative phase in the DE-neutrino 
coupling matrix, respectively. There is also the usual 
vacuum mass difference (squared) /S.m? and the vacuum 
mixing angle 9^. The analogous propagation equations 
for antineutrinos are the same as Eq. ([5]) with the re- 
placements moff — > — TOoff and rj — > —77, the changes in 
sign a manifestation of CPT violation. 

Recall that if the propagation Hamiltonian is of the 
form, 

/, = Mf-^^°!f '"^Z], (6) 
\ sm2t^ cos26' J ' 

I 



then the probability for one species of neutrino to convert 
to a different neutrino after a distance L is 



Pii^a -> i^b) = sin^ 29sm^{ML). (7) 



Here we have neglected the CP- violating phase in Eq. ([6]) 
because it does not affect the oscillation probability. The 
propagation Hamiltonian in Eq. ^ can be written in the 
form of Eq. ^ with the following relations [22] : 



sin^ 261 



4:E 



mcs — V ■ p) Am^ 



AE 



"^off (1 — V ■ p)(cos 29y cos 29d + sin 20^, sin 29d cos 77), 



1 

M2 



4:E 



sin-^ 29„ 



mcff^(l - t> -p) . 



sin-" 29d 
I 



Am 

iE 



TOoff (1 — V ■ p) sin 26*1, sin 29 d cos rj 



(8) 
(9) 



Note that this time sin 29 does indeed depend on 77, as it not be rotated away by redefinition of wave functions, 
is not the overall phase, but the relative one, that can- 
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The oscillation length is then Lose = ttM~^. In the ab- 
sence of DE-neutrino coupling, we recover the standard 
oscillation length Lose = 47ri?/Am^ and mixing angle 
6 = 9y. If nicB ^ Am^/2E, then the oscillation length 



is 



Lose = 27rTOcff (1 -v-pY 



In general, rrics can be either positive or negative. 
However, from the symmetry of the Hamiltonian, the 
relevant parameter space can be limited to rrice > 0, 
<9d< 7r/4, and < 77 < 77 

Thus far, no deviations from standard three- flavor neu- 
trino oscillations have been discovered in experimental 
data (except LSND [lO]), and this yields constraints on 
CPT-violating parameters, especially for meff- By ana- 
lyzing solar-neutrino and KamLAND data, Ref. ob- 
tained an upper limit of nice < 3.1 x 10"^'' GeV for i^g-t'^ 
mixing. Atmospheric and accelerator data provide an up- 
per limit for Vfj^-Vr mixing of rries < 5 x 10~^^ GeV |14| . 



III. ULTRA-HIGH-ENERGY NEUTRINOS 

A. Prediction 

Given that DE-induced neutrino mixing becomes in- 
creasingly important, relative to vacuum mixing, at high 
energies, the DE-neutrino coupling can be probed with 
ultra-high-energy cosmogenic neutrinos. These neutri- 
nos are produced by the interaction of ultra-high-energy 
cosmic-ray protons with CMB photons [1^ : 



PI - 



(10) 



The fact that the Greisen-Zatsepin-Kuzmin cutoff 24 1 
has now been observed by the HiRes [IHl and Auger 26 j 
collaborations implies that this interaction must be oc- 
curring. And if so, there must be a population of cosmo- 
genic neutrinos with energies 10^^-10^° eV [1^, [H, 

The characteristic distance between the source of these 
neutrinos and the Earth is the Hubble distance cH^^ ^ 
which is much longer than the oscillation length — i.e., 
cHq'^ > A/"i~as long as rrics > iJo = lO"''^ GeV, as 
is always the case here. Therefore, any oscillatory fea- 
tures in neutrino mixing will be washed out; the proba- 
bility for conversion of a cosmogenic neutrino from its 
production flavor to another flavor en route from the 
source is then simply sin^26'/2. Cosmogenic neutrinos 
mostly originate from pion decays, with the characteris- 
tic flavor ratio v^w^Wt — 1:2:0. The result of standard 
vacuum mixing would be a flavor ratio at the Earth of 
Ve-v^-.v.,- = 1:1:1. While possible corrections to this fla- 
vor ratio can be induced by small three-flavour oscillation 
effects or other new physics, here we concentrate on ex- 
ploring the consequences of the DE-induced mixing. 

For the sake of simplicity, we focus on v^-Vr mixing 
(and their antiparticles) . In the absence of a DE-neutrino 
interaction, these two flavors are maximally mixed; i.e., 
Ov = 7''/4, and thus even if only are produced at the 



source, an equal number of and Vt is generated by 
mixing. However, this can be altered if there is a DE- 
neutrino interaction. The flux of and v^- at the detec- 
tor is related to the flux at the source through. 



= I 1 - i sin^ 29 



■sin^ 29 . 



(11) 
(12) 



and will in general differ from 9^ if meff 7^ 0. 



B. Proposed Measurement 

Here we investigate the possibility of measuring 9 using 
current or future ultra-high-energy-neutrino experiments 
such as Auger 30] and ANITA [3l| . It is in principle pos- 
sible to discriminate Vr from j/^ by separately measuring 
the Earth-skimming events (^'t) and almost horizontal 
events originating in air (v^). 

We assume that a given experiment detects iV^°* neu- 
trino events. This quantity is for the total neutrino and 
antineutrino flux; i.e., Njf^ — + N^^ (here v rep- 
resents both neutrinos and antineutrinos). If the fla- 
vor democracy expected from vacuum mixing is realized, 
then one expects Ny^ = Ny^ = A^*°*/2. The number 
of neutrino events at the detector is related to the flux 
through. 



Ny 



(13) 



where S(£^) is the detector exposure to neutrinos in units 
of cm^ s sr, and it generally depends on neutrino energy. 
Here we assume 0*°* = (p^ = KE^^ with a normaliza- 
tion constant K, L^min = 2 x 10^'' eV, and i^max = 2 x 10^^ 
eV. This provides a good approximation for the spectrum 
of cosmogenic neutrinos (e.g., Ref. [29]). For simplicity 
we further assume that the detector exposure is indepen- 
dent of energy; the Auger ex pos ure indeed depends on 
neutrino energy only weakly [30|. Therefore, the total 
number of neutrino events is given by 



KE 

-£■111 in 



(14) 



and the number of Vr events is given by 



Ny 



dE -sir? 29(1)1 (E)^ 



KE 



E„ 



dE E-^ sin^ 261 



tot TP pEj 



NTE, 



dEE-^sm'^29, (15) 



where we used Eq. (fHl) in the last equality. 
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FIG. 1: Sensitivity on (mcff, sin^ 2S(j) plane of future exper- 
iments that would yield Nl°^ = 12 and 100 total neutrino 
events. 



To investigate the sensitivity of a given experiment, we 
assume a null detection of new physics; i.e., the result of 
Ni,^ is consistent with the standard expectation N^"^ /2 
within statistical errors (we do not take systematic un- 
certainties into account). This will reject a certain range 
of parameter space for {irieS, sin^ 29d). More specifically, 
to obtain 95% C.L. (2cr) limits for these parameters, we 
solve 



(16) 



for nics and 6, using Eq. (fT5)) for the left-hand side. In 
Fig. [1] we show the sensitivity of detectors that are ex- 
pected to collect 12 and 100 neutrino events^ (total) and 
that also have a j/T--identification capability. If the true 
values of rrics and sin^ 29d are above these curves, then 
we will see an anomalously suppressed Vr flux compared 
with the standard mixing scenario. We also show the 
current upper limit on rrios obtained from the combined 
analysis of Super-K and K2K data performed in Ref. [l3| . 
One can see from this Figure that by detecting cosmo- 
genic neutrinos and by studying their flavor content, one 
can largely improve the current sensitivity to m^s and 
9d, quantifying further the suggestion of Ref. [l3|. We 



^ The current Auger exposure is ~ 10^'^ cm^ s sr |3Cli] . and an 
optimistic estimate for the fiux of cosmogenic neutrinos is close to 
the Waxman-Bahcall bound [H], E^(j>^{E) ~ 10"* GeV cm'^ 
s-l sr-l m. Therefore, from Eq. ((HJ, we expect A^^"' < 1, 
which is still consistent with nondetection by Auger. 



also note that a weaker sensitivity, albeit still much bet- 
ter than the current sensitivity, may be achieved with 
neutrinos of slightly lower energies [32| . 



IV. MATTER EFFECTS IN ATMOSPHERIC 
NEUTRINO OSCILLATIONS IN THE PRESENCE 
OF A DARK-ENERGY COUPLING 

We now turn our attention to Lorentz/CfT- violating 
effects in electron-neutrino oscillations, showing here that 
novel effects may arise with DE-neutrino coupling as neu- 
trinos propagate through the Earth. These effects may 
allow us to access with atmospheric neutrinos regions of 
the DE-ncutrino-coupling parameter space significantly 
below those currently probed. In this section, we consider 
two-flavor and three-flavor oscillations. 

As neutrinos travel through matter, there is an ad- 
ditional contribution to oscillations from the matter po- 
tential \/2GFNe (where Gf and Ne are, respectively, the 
Fermi constant and electron density) relevant if electron 
neutrinos are involved. Recalling that the matter poten- 
tial is > 10~^^ GeV, the vacuum-mixing term /^m? /2E 
is small for neutrino energies > 10 GeV. The mixing ma- 
trix Eq. ([5]) then becomes for v^-v^ mixing (neglecting 
the overall factor of 1/2, the directional dependence, and 
the phase rf), 



—rrics cos 

meff sin 29 d 



rrics sin 29d 
rrics cos 29d - V2GfN, 



(17) 

Note that here, both the DE term and the matter poten- 
tial change sign for antineutrinos, unlike the usual MSW 
effect, in which the vacuum term does not change sign. 
Unlike MSW mixing, there is essentially no energy depen- 
dence, at sufficiently high energies, in this mixing matrix. 

To see when DE- induced mixing may be significant, re- 
call that the value of the matter potential is \/2GFNe — 
7.6 X 10^^'^Ye{p/g cni~^) eV. The Earth core has av- 
erage density pcorc = 11.83 g cm~^ and electron frac- 
tion Y^°'^'^ = 0.466, while the mantle has average den- 
sity Pmantic = 4.66 g cm'^ and y/i'^itie = 0.494, with 
the surface layer of the Earth having density as low as 
2.6 g cm^"^. The matter potential is thus about 10^^'^ 
eV, so the effects of DE-induced mixing may be manifest 
for TOcff around 10^^^ GeV, well below current upper lim- 
its. In the absence of matter, as discussed in Ref. [ll|, it 
is possible to obtain a resonance when all mixing angles 
involved are maximal 



Am 
2E 



cos 29v + cos 29d = 0. 



(18) 



Here, in the presence of matter and at high energies, a 
resonance can occur for a small mixing angle 9d when 



cos26'rf = V2GFNe. 



(19) 



The presence of a resonance is thus entirely determined 
by the densities encountered along the path and the DE 
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FIG. 2: Oscillation probability as a function of zenith angle 
for atmospheric neutrinos of i5 = 50 GeV, obtained with 9d — 
n/4. 



rrics = 5 X 10"^"^ GeV for three-flavor mixing in the 
Lorentz-violating sector. The case where = 7r/4 and 
= also corresponds to an effective two-flavor sce- 
nario, just like the previous results. It leads, however, 
to a very different behavior due to the different contribu- 
tion of the standard neutrino oscillations. The case where 
^13 — ^12 = 7i'/4 corresponds to a full three-flavor oscilla- 
tion scenario. We have also studied the effects for other 
values of the mixing angles and the same features remain 
present. A non-zero value of for standard neutrino 
oscillations leads to similar features in the zenith-angle 
distribution. However, the effects are extremely small 
at the high energies considered here, orders of magni- 
tude below those coming from the Lorentz-invariance- 
violating terms. 



DIRECTIONAL DEPENDENCE 




FIG. 3; Oscillation probability as a function of zenith an- 
gle for atmospheric neutrinos of iJ = 50 GeV, obtained for 
three-flavor oscillations 
mcff = 5 X 10"^^ GeV. 



coupling parameters, with no (or very weak) energy de- 
pendence at high energies. 

To illustrate the possibilities, we integrate the 
neutrino-propagation equation (including the small 
vacuum-mixing term) to calculate the Vfj^-to-Ve transi- 
tion probability as a function of (cosine of) the zenith 
angle for atmospheric neutrinos propagating through the 
Earth. We use the de nsity proflle of the Earth as given 
by the PREM model (S^l • Figure [H shows the results 
for two-flavor oscillations for different values of mcfr for 
0d = tt/4:. When rrics ^ V^GpNe, the oscillation prob- 
ability is determined almost entirely by the DE term; 
there are regular large-amplitude variations of the oscil- 
lation probability as a function of zenith angle. As rrics 
decreases to values comparable to v^G^iVe, the oscilla- 
tion probability decreases, and the oscillation length is 
seen to differ for trajectories that do {cos 9 < —0.8) and 
do not (cos 6* > —0.8) pass through the core. 

In Fig. [3] we show the oscillation probabilities for 



While detection of CPT/Lorentz-violating effects 
would be spectacular — it would imply new physics re- 
gardless of whether it is DE-related or not — the real 
"smoking gun" for a DE effect would be the directional 
dependence, oc {1 — v -p), of neutrino-oscillation param- 
eters. Given that our peculiar velocity with respect to 
the CMB rest frame is 10^"^ times the speed of light, the 
magnitude of this effect is going to be suppressed relative 
to the other effects, discussed above, of a DE-neutrino 
interaction. Statistics well beyond the reach of current 
and forthcoming neutrino experiments will be required 
to detect this effect. Still, it is worth keeping in mind for 
future generations of experiments. 

It may also be worth searching for such a directional 
dependence in current data, just in case there is a DE- 
neutrino coupling that is manifest in ways different than 
we have foreseen here. For example, if DE somehow pro- 
duces Lorentz violation through a modiflcation of the ki- 
netic term in the Dirac equation, the energy dependence 
of the mixing induced by Lorentz/ CPT- violation could 
be different We therefore work out in this Section 

expressions for the factor v ■ p to aid experimentalists 
who may wish to look for direction-dependent effects in 
their neutrino (or other) data. 

To proceed, we first set our coordinate system. We set 
the origin at the center of Earth and align the z axis along 
the rotational axis of Earth, so that the north pole has 
positive z coordinate. We set the x axis along the direc- 
tion to the Sun at vernal equinox. Since the Sun moves 
eastbound, its position at summer solstice aligns with the 
y-axis. We can thus represent the seasonal shift by an az- 
imuthal angle tp, where ip — 0, tt/2, tt, and 37r/4 for vernal 
equinox, summer solstice, autumn equinox, and winter 
solstice, respectively. Note also that the orbital plane of 
the Sun is inclined from the x-y plane by 9inc — 23.5° [34| . 

The Sun is moving with respect to the CMB rest 
frame with a speed of vq = 369 km s~^ towards 
the direction a = 168°, 6 = -7.22° where a 

is right ascension and S is declination of the celes- 
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tial coordinates [3J|. In our coordinates, the velocity 
of the Sun is Vq = Wq (cos 5 cos a, cos 5 sin a, sin J) — 
(—358, 76.1, —46.4) km s^^. The Earth is moving around 
the Sun with average orbital speed oi V = 29.8 km s~^. 
Thus, the velocity of the Earth with respect to the CMB 
rest frame is 



-358 

76.1 
-46.4 



simp 

— cos ip cos Oinc 

— COS if sin 9inc 

+29.8 sin 
— 27. 3 cos 

— 11.9 cos(p 



km s 



(20) 



We neglect the contribution from the rotation of Earth 
(< 0.5 km sec~^) to our velocity with respect to the CMB 
rest frame. 

Now we evaluate the direction of the neutrino beam 
p. We suppose that the beam runs from some point A 
on the Earth's surface to another point B on its surface 
(or vice versa). It should be straightforward to general- 
ize the arguments below so that extraterrestrial neutrino 
production can be taken into account. We set the ori- 



gin of time coordinate T to "noon" (i.e., when the Sun 
reaches highest) at the point A. Therefore, the positions 
of A and B in our coordinate are 



cos (j)A cos(a;rA + f) 
xa ~ Rq I cos (j)A sin(u;TA -I- if) 
sin (pA 

cos 4>B cos{lljTa + AA 
Xb = I cos (j)B sin(a)TA -f- AA 
sin <j)B 



(21) 



(22) 



where i?0 is the radius of the Earth, w is the rotational 
frequency (27r/day), Ta is the time at the position A 
relative to noon, 0a, b is the geometric latitude of the 
points A and B, and A A = Aa — Ab is the difference of 
the geometric longitude. The quantity AA appears in xb 
because we measure the time (for both A and B) with 
respect to noon of the point A, so the time difference 
is given by the latitude difference (note also that the 
longitude increases to the west). The direction of the 
neutrino beam p is then proportional to xb — xa with 
proper normalization as 



p = 



■\/2(l — cos (fiA COS (f)B COS AA — sin (f)A sin (f)B) 



cos <j)B cos{ujTa + AA + f ) — 
COS <j)B sin(Li;TA + AA 



- AA + If) — COS 0A cos{ujTa + v)\ 

- AA + (f ) — COS (j)A sin(a;TA -t- f) 
sin ~ sin <^a / 



(23) 



Therefore, by combining Eqs. (|20p and ([23]), we obtain 
the directional factor 1 — ■ p. Since it is a scalar 
quantity, the final result does not depend on the choice 
of the coordinate system. 



VI. THEORETICAL IMPLICATIONS 

Before closing, we discuss, for illustration, the im- 
plications of a measurement of a particular value of 
the Lorentz-invariance-violating effective mass parame- 
ter rrieff in terms of a specific model of DE-neutrino cou- 
pling. 

Perhaps the simplest interaction of this kind has the 
form, 

Cint = -Aa/3-^j>a7''(l " 15^13: (24) 

where is a quintessence field, Xap is a coupling- 
constant matrix, and Af* is some mass scale. Thus, 
a1 ~ A^(t)Z''/M*, and rjias ~ AA0(i)/M*, where AA is 
the difference between eigenvalues of the A matrix. For 
quintessence, one expects 4' ~ MpiHo{l + w)^/^ (e.g., 
Ref. 3), where Mpi is the Planck energy scale. In this 



case, the mass scale M.^, corresponding to a given TOoff is 

to the mass scale that controls the DE-neutrino in- 
teraction. The ultra-high-energy Vf^r^r oscillation ef- 
fects we have discussed thus probe up to mass scales 
Af* ~ 10^ GeV. The I'e-t'^ oscillations induced by the 
matter effects we discussed probe up to mass scales 
~ 100 MeV. 



VII. CONCLUSIONS 

We studied the implications of an interaction between 
dark energy and neutrinos for neutrino oscillations. The 
most general Lorentz/CPT- violating term induced by 
dark energy (DE) takes the form {aLYv^^{l — 75)1^, 
where [aL)^ is a four- vector normal to the CMB rest 
frame. This introduces a new source for neutrino os- 
cillations that are energy independent and different for 
neutrinos and antineutrinos. Furthermore, the motion of 
the Earth with respect to the cosmic rest frame induces 
a directional dependence in the oscillation probabilities. 
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The current best limits to the DE-neutrino cou- 
phng we considered are obtained from atmospheric- and 
accelerator-neutrino experiments for Vf^-Vr mixing, and 
from solar and reactor experiments for Vf^-v^j_ mixing. 
However, the higher the neutrino energy, the more promi- 
nent the effect of the DE-neutrino interaction. We 
therefore considered in this paper cosmogenic ultra-high- 
energy (energies of 10^^-10^^ eV) neutrinos produced 
by the interaction of ultra-high-energy cosmic rays with 
CMB photons. We showed that future experiments tar- 
geting these neutrinos will improve the sensitivity to a 
DE-neutrino interaction by seven orders of magnitude, 
down to mcff ~ 10^^° GeV compared with the current 
upper bound mcff < 5 x 10^^^ GeV (Fig. [1]). This corre- 
sponds to a sensitivity to an energy scale as large as ~10^ 
GeV for the DE-neutrino interaction. We then showed 
that the interplay of DE- and matter-induced neutrino 
mixing could induce a novel zenith-angle dependence for 
i^e oscillations in atmospheric neutrinos. This effect may 
extend the sensitivity to Lorentz/CJ-'T'- violating param- 
eters in the v^. by roughly three orders of magnitude. 

The real smoking gun of a DE-neutrino interaction (as 
opposed to some other origin for Lorentz/CPT violation) 
would be a directional dependence of the oscillation prob- 
abilities. The notion that Lorentz violation may give rise 
to a directional dependence is not new (e.g., Ref. [36|) 
and searches for directional dependence in neutrino ex- 
periments have already been carried out (e.g., Ref. [13]), 
but prior work has considered Lorentz-violating param- 
eters introduced in an ad hoc manner and/or tested 
for direction-dependent effects in a Sun-centered inertial 
frame. We emphasize here that cosmic acceleration sug- 



gests that we seek a specific form of Lorentz violation, 
that where the preferred frame is aligned with the CMB 
rest frame. Even though such a signal is expected to 
be small, it is still worth seeking in existing and future 
experimental data. 

We have not discussed specific models for a DE- 
neutrino interaction, beyond an illustrative toy model, 
but it may be interesting to do so (see also Ref. [ssj). 
The theoretical motivation to expect such a coupling may 
admittedly be slim. However, we are at square one in 
our understanding of DE, and such a coupling is no less 
likely to be expected, perhaps, than any of the many 
other manifestations of new cosmic-acceleration physics 
that have been considered. Discovery of Lorentz/ CP T- 
violating effects would be extremely important, even if 
not attributable directly to dark energy. A directional 
dependence, if discovered, would be absolutely remark- 
able, as it would provide moreover clear evidence that 
there is more to cosmic acceleration than simply a cos- 
mological constant. 
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